Time Delays in the Synchronization 
of Chaotic Coupled Systems with Feedback 
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Isochrony and time leadership was studied in the synchronized excitable behavior of coupled 
chaotic lasers. Each unit of the system had chaos due to feedback with a hxed delay time. The inter- 
units coupling signal had a second, independent, characteristic time. Synchronized excitable spikes 
present isochronous, time leading or time lagging behavior whose stability is shown to depend on a 
simple relation between the feedback and the coupling times. Experiments on the synchronized low 
frequency fluctuations of two optically coupled semiconductor lasers and numerical calculations with 
coupled laser equations verify the predicted stability conditions for synchronization. Synchronism 
with intermittent time leadership exchange was also observed and characterized. 
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Many dynamical systems in nature have chaos due to 
feedback. A characteristic time, tf, associated to the 
feedback is therefore imbedded in the system response. 
As two or more of such systems are coupled, another in- 
dependent time, Tc, corresponding to the time of flight of 
the coupling signal, enters in the dynamical description of 
the global system. We present here how the relation be- 
tween these two times determines the possible time delays 
in chaos synchronization. New features of time leader- 
ship competition in synchronized chaos between coupled 
pairs of systems with feedback are found when feedback 
and inter-coupling times have the same order of mag- 
nitude. Experimental and numerically, the stability of 
isochronous chaos synchronism for identical systems oc- 
curs only for specific relations between these times. We 
also show that time delayed and time advanced synchro- 
nism as well as synchronism with intermittent leadership 
exchange are also quantitatively determined by the ratio 
of these times. Our case is made with pairs of semicon- 
ductor diode lasers. However, the properties of synchro- 
nization in complex systems extends far beyond physical 
devices reaching the subject of neural sciences Q. 

The study of synchronism with chaotic lasers spreads 
for more than a decade [1, 0) S S 0] ■ Of relevant interest 
for applications are the results on the synchronization of 
semiconductor lasers [1,0] for encrypted communication. 
With optical feedback diode lasers can present chaos in 
the form of very fast output power fluctuations, at a time 
scale of picoseconds. Superimposed on these, irregular 
power drops occur in a much slower time scale (order 
of hundreds of nanoseconds and longer) corresponding 
to the so called low frequency fluctuations (LFF) [lo| . 
The LFF drops in single lasers are a current subject of 
studies and have been associated to spikes of excitable 
systems [lli, Coupled diode lasers show time 

advanced and time lagging synchronization via unidirec- 
tional coupling in master-slave configuration [3, and 
in mutually coupled systems [13] . Symmetrically coupled 
pairs of lasers, without feedback, were shown to have un- 



stable isochronous chaotic pulsation [3, El, El. In the 
experiments and calculations with coupled lasers without 
feedback, the time leading lasers always appears with its 
power drops displaced by one unit of the coupling time, 
Tc ■ The use of intermediate relaying system was demon- 
strated to give stable isochrony |l^, l20j. Isochronous 
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synchronization has also been investigated [2C 
for lasers with feedback with the studies focused on the 
fast laser fluctuations. Differently, herein we study the 
synchronization of the low frequency fluctuations (LFF), 
known to appear at the scale of hundreds of nanoseconds 
and slower. Thus, in our case the dynamics in the cou- 
pled systems is much faster than the synchronized events 
which have the properties of excitable spikes. So, our 
results refer to synchronism of excitable dynamical sys- 
tems. The schema of the experiments is given in figure [T] 
Optical feedback was created in each laser by a retrore- 
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FIG. 1: Experimental setup for the synchronization of Low 
Frequency Fluctuations in optically coupled laser 

fleeting external mirror. Their feedback return times, 
Tpi and Tp2, were set equal to within 1% precision (both 
named from hereon Tp). The time of travel of the cou- 
pling signal between the lasers, tq, was independently 
controlled with respect to tf- Small changes in either of 
the times, on the scale of fraction of nanoseconds do not 
alter the properties of the synchronism. Thus the results, 
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like the LFF phenomenon in single laser with feedback, 
are robust with respect to optical phase changes. A con- 
sequence of such behavior is that our observations are 
consistent with the synchronization of excitable systems 
[2^. These authors show how unidirectionally coupled 
excitable systems can presents one of the systems always 
with a lower threshold for excitability. Their coupled sys- 
tems synchronize to an external common signal. In our 
case the coupled systems have fast fluctuations and so, 
no external source of excitation is necessary. 

The experiments were done with two SDL 5401 
GaAlAs semiconductor lasers, named here as Laser 1 and 
Laser 2, both with solitary threshold current of 21 mA 
and emitting at 805 nm. They were thermally stabilized 
to 0.01 K and could be temperature tuned to have their 
optical frequencies within 2 Ghz separation. The feed- 
back time of the lasers was tf = 10 ns and the strength 
were measured by the threshold reduction parameter, ^ 
which is the percentual variation of the threshold pump 
current as we consider the laser with and without feed- 
back [1^. Symmetrical optical coupling between the 
lasers was produced with 50% beam splitters as shown in 
figure [H The time of flight of the light between the lasers 
was varied between 5 and 20 ns. The threshold reduc- 
tions due to cross input power were used to quantify the 
coupling strengths whose contribution to our studies will 
be detailed elsewhere. Manipulating the laser currents 
and temperatures, LFF synchronism was obtained where 
to each power drop in laser 1 corresponds a drop in laser 
2 and vice versa. Typical experimental segments of the 
power output of the two lasers, with three events of the 
synchronized low frequency fluctuations (LFF) drops, are 
shown in figured Each laser intensity was detected by a 
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FIG. 2: Segment of experimental time series of the power of 
(a) laser 1 and (b) laser 2, showing the low frequency synchro- 
nism and including the interchange of time delay between the 
lasers, tf = 10 ns and tc = 10 ns. 

1.5 GHz bandwidth photodiodc. Simultaneous data se- 
ries were acquired with a two channel digital oscilloscope 
having a bandwidth of 1 GHz and a maximum sampling 
rate of 5 GSamples/second. One of the lasers could al- 
ways be the time leader by setting its pump current or 
its optical frequency [l6| higher. However, as careful tun- 



ing of the laser currents and optical frequencies is made, 
within the same long LFF synchronized time series the 
LFF drops of laser 1 can appear isochronous, leading 
or lagging these of laser 2. Such were the conditions of 
figure [5] where the time leadership is interchanged be- 
tween drops one and three, while the second LFF drops 
are isochronous. Using long pairs of data series, coarse 
grained to a 1 ns time resolution [l3| . we measured ATi— 
Tii — T2i, the delay between the LFF drop i of laser 1 
and laser 2. 

The main result of this work is to show that, within 
any LFF synchronized evolution, for any pair of excited 
drops the allowed delays are given by 



AT!; = rui ■ Tc - Ui ■ Tp 



(1) 



where ^ and rii are small integers. Here we only 
give indications of m,; = ±1 but preliminary results with 
very long numerical series show rare events with > 1. 
Equation [1] covers various previously studied cases in the 
literature and confirms novel observations. For instance, 
with the lasers having no feedback, tf = 0, it verified 
that isochrony is not allowed [3] ■ Other cases with feed- 
back but forbidding isochrony are given below. The equa- 
tion determines ATi when one laser is always leading or 
lagging which happens for a single pair {mi^rii) through 
the whole time series. Furthermore, it is also valid for 
cases with varying values of (m^, Ui) within the same se- 
ries, obtained from symmetrical and nearly symmetrical 
systems. The LFF synchronism then appears with inter- 
mittent switching of the time delay that can interchange 
the leading subsystem. The evidence of equation [1] as a 
property of the LFF synchronism was obtained inspect- 
ing many different feedback and coupling delay times in 
experimental and numerical series. The origin of equa- 
tion [1] is present in the intensity cross-correlation func- 
tions of the fast fluctuations of coupled lasers. These 
correlations, obtained with sub-nanosecond resolution, 
show recurrent narrow spikes (width of hundred of pi- 
coseconds) at positions and intervals given by[TJ 

Let us proceed presenting our experiments along with 
the numerical-theoretical results extracted from a sys- 
tem of differential equations that describe two mono- 
mode lasers with optical feedback and optical coupling. 
The model correspond to a set of modified Lang and 
Kobayashi [23| equations for the lasers with feedback, 
including symmetrical optical coupling and assumed to 
have the same optical frequency: 
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where i,j = 1, 2 and each laser gain is given by 

GoiN,-No) 



1 



mt)\ 



(4) 



The definition of the various parameters and their typ - 
ical values are well discussed in the literature [l^, 
For each laser, Ej is the radiation field, a is the factor 
describing amplitude to phase conversion, G is the ampli- 
fying gain, Nj (t) the carriers inversion population, Tp the 
photon lifetime of the internal laser cavity, Tg the carriers 
lifetime, Jj the threshold normalized pump currents and 
A^o the inversion population for medium transparency. 
Each feedback field has an amplitude coefficient k and 
feedback time Tp. The optical couplings are linearly ad- 
ditive E field with coefficient 7 and time tq for the field 
of one laser to reach the other one. Physical causality 
demands that both Tp and tc be positive. 

With a fourth order Rungc-Kutta algorithm numeri- 
cal data series were obtained for Ei (t) and these gave the 
normalized intensity series \Ei{t)\'^. The time scales in 
the integrations and the equations parameters were fixed 



as 

K= 



dt=l ps, l/rp=282 ns-\ A^o 
=7 = 22 ns-\ 1/ts = 1.66 ns~ 



1.5xl0«, e = 5x10" 



The times, Tp and tc, 
are in the nanosecond range. Calculation with one laser 
having sufficiently higher pump current gives synchro- 
nized LFF with the laser of higher pump current always 
leading in time, as observed in the experiments. Equal 
pump currents (Ji = 1.013), corresponding to symmetri- 
cal systems, produced numerical time series with leader- 
ship exchange, again in agreement with the experiment. 
Segments of these data series are shown in figure [H to 
be compared with the experimental segments in figure [2l 
Let us emphasize that no stochastic terms are used in 
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FIG. 3: Typical segment of tlie numerical time series of (a) 
laser 1 and (b) laser 2 LFF power drops under synchro- 
nized condition. Notice that the first drops are isochronous 
while time leadership was interchanged between drops two 
and three, tf ~ 10 ns and tc ~ 10.05 ns. 



the equations. The apparent random time distribution 
of the LFF events in each laser as well as the time lead- 
ership switching in the synchronism is due to the high 
dimension deterministic chaos. 



Time series with ~ 10*^ pairs of drops were coarse grain 
filtered [13| with a 1 ns time constant, and used to ex- 
tract histograms of time delays between the drops. The 
switching of leadership from symmetrical systems was 
also examined and no correlation was found between con- 
secutive values of AT^, up to second order conditional 
probability. This is indication of a Markov process, ob- 
tained despite the fact that the data came from determin- 
istic numerical equations with recurrence times Tp and 
Tc ■ Such result is consistent with high dimensional chaos 
with largely different time scales (Ti+i — Ti ^ tp,tc) 
that makes LFF appear as stochastic without memory 
[ist . The experimental data show the same lack of cor- 
relation for the delay times of LFF synchronism in sym- 
metrical conditions. The histograms associated to the 
data series of figures [D and [31 are shown in figure |31 As 
Tp = 10 ns and tc — 10 ns, equation[T]with rm = ±1 and 
rii = ±1, predicts ATi — 0, ±10ns and ±20 ns. Indeed, 
the histograms show that isochrony occurs, along with 
time leadership exchange events. The major probabili- 
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FIG. 4: Histograms of the (a) experimental and (b) numeri- 
cal time delays between LFF pulses of the two synchronized 
lasers, tf = 10 ns and tc ~ 10 ns. 

ties are for events with AT; — and ± 10 ns, with few 
events at ±20 ns and some in ±30 on the experimental 
data. It is important to mention that on the calculations 
small differences (less than 2%) were introduced between 
the parameters of laser 1 and 2 in the equations. In par- 
ticular Tc ^ Tp [tc — 10.05 ns) was necessary to give the 
non isochronous events. In fact such differences which are 
intrinsically present in the experiments change dramat- 
ically the amplitude of the peaks in the histograms but 
have minor effect on the values of allowed delays. The 
robustness of the relation between delay times and the 
condition of LFF synchronization according to equation 
[1] was also inspected for small, up to 10%, of relative 
variations of Tp and tc- 

Another typical numerical histogram where rui and rii 
change within the same synchronized evolution is shown 
in figure \5\ This is also a case of symmetrical system 
where the laser parameters are equal but the values for 
time of coupling and feedback, were chosen to give unsta- 
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FIG. 5: Histograms from experimental (a) and numerical (b) 
series for the time delays between LFF drops of the two syn- 
chronized lasers: tf = 10 ns and tc = 15 ns. 



ble isochrony: tc = 15 ns and Tp = 10 ns. Accordingly, 
equation [1] these values prevent AT^pp — 0, as seen by 
direct substitution of small integers (m^ — ±1). Thus, 
synchronism of the LFF holds, but there is always a fi- 
nite time lead between the two lasers. The leadership 
can be exchange but events of simultaneous drops are 
excluded. 

When Tc ^ Tp we verified numerically that equation 
[T] still holds. The value of is always ±1 while rii 
assumes a large range of values. Isochrony is absent and 
the dominant events occur with rii = 0, corresponding to 
delays of ±tc- Nonsymmetrical systems, as mentioned 
above, also follow equation [TJ A calculation with the 
two lasers having the same current Ji = 1.013 was given 
figure [S] (a) while figure [S] (b) shows what happens when 
Laser 1 has its current augmented to Ji = 1.025. The 
laser with higher current is always the time leader, even 
though sometimes its leading time jumps down from tc 
to Tc — Tp. These events correspond to be = 1 and 
Hi — changing to = 1 and rii — +1 along the 
dynamical evolution with synchronized LFF excitations. 
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FIG. 6: Histograms of the numerical delay times between 
synchronized LFF drops of the lasers with tf = 10 ns and 
TC ~ 15 ns and showing the effect of asymmetric conditions. 
In (a). Ji = 1.013 for the two lasers and in (b) Ji — 1.025 
and J2 = 1.013, making Laser 1 the time leader. 



To summarize, from experiments with coupled lasers 
with feedback, corroborated by the numerical solutions 
of the corresponding rate equations, we discovered that 
the delay time in the synchronization of coupled ex- 
citable systems with feedback is controlled by a sim- 
ple relation between feedback time and inter-coupling 
time. One of the systems may have a fixed time lead- 
ership or have its leading time switching values including 
isochrony. Nearly symmetrical systems can also have in- 
termittent time leadership exchange always maintaining 
the excitable spikes synchronized in the large time scale. 
Specifically, the behavior of identical systems with feed- 
back and coupled with time delay strongly depends on 
the relation between these times. A simple equation that 
specify the intercombinations of feedback and coupling 
time to give the allowed values of delay in the synchro- 
nized dynamics was introduced. Coupled asymmetrical 
systems also follow the conditions for the allowed delay 
time between events. The various parameters of the sys- 
tems strongly influence the probability of specific delays, 
distorting the histogram amplitudes, but have no effect 
on the values of the time delays. Our results extend 
the previous [l^ determination of the symmetry break- 
ing and instability of expected isochrony in the synchro- 
nism of identical coupled chaotic systems without feed- 
back. It also adds to the recently studied properties of 
isochronal synchronism done by Klein et al. 21| and Kan- 
ter et al. 0], who emphasize the potential application of 
these properties for encrypted communication. A formal 
mathematical treatment of the stable and unstable man- 
ifolds in chaotic synchronization with delays having in- 
teger combinations between feedback and coupling time 
will be presented elsewhere. The phenomenon of selec- 
tion condition in time delays is bound to appear generi- 
cally in mutually coupled dynamical systems. 
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